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Control of dispersion in fibre optical waveguides 
is of critical importance to optical fibre commu¬ 
nications system^^ and more recently for con¬ 
tinuum generation from the ultraviolet to the 
mid-infrarecPM^. The wavelength at which the 
group velocity dispersion crosses zero can be set 
by varying fibre core diameter or index stepP^^®^. 
Moreover, sophisticated methods to manipulate 
higher-order dispersion so as shape and even flat¬ 
ten dispersion over wide bandwidths are possible 
using multi-cladding Here we introduce 

design and fabrication techniques that allow anal¬ 
ogous dispersion control in chip-integrated op¬ 
tical microresonators, and thereby demonstrate 
higher-order, wide-bandwidth dispersion control 
over an octave of spectrum. Importantly, the fab¬ 
rication method we employ for dispersion control 
simultaneously permits optical Q factors above 
100 million, which is critical for efficient opera¬ 
tion of nonlinear optical oscillators. Dispersion 
control in high Q systems has taken on greater 
importance in recent years with increased interest 
in chip-integrable optical frequency comb^I^H^. 


High-Q microresonators, through cavity enhancement 
of an input field, enable many important optical devices 
and functionJ^. These applications i nclud e frequency 
comb generatioiP^M^, pulse source^^^^, casc aded 
Raman laserJ^^i^, stimu lated Brillouin laserJ^HlMl and 
harmonic generatioiP^^^. In many systems the control 
of modal dispersion is helpful or essential. This is par¬ 
ticularly important in parametric oscillator d^^ l ^^ l where 
the local dispersion must be anomalous and in frequency 
microcombs where both the sign of the dispersion a s well 
as its spectral shape are critical for comb operatiorP^^^. 
Spectral bandwidth and coherent operation are strongly 
influenced by dispersion, and there has been remark¬ 
able progress in these areas, including sub-co mb 
synchronizatiorP ^ I ^^ I ^^ I ^^ I, so liton gene ratior l^^ l ^^ l ^^H ^, 
and dispersive wave formatiorP^^^^^^^^. 


There have been multiple approaches to engineer 
dispersion in microresonators. The control of waveguide 
width and height in 8 ^X 4 , diamond, MgF 2 , an d CaF 2 
resonators alters geometric dispersioiP^l^^^^^l^. Hf 02 - 
coated Si 3 N 4 an d oxid ized Si resonators also provides 
dispersion controP^^. This paper studies dispersion 
control in ultra-high-Q (UHQ) silica disks using a 
method that is inspired by dispersion engineering in op¬ 
tical fibres. By creating multi-wedge disks and precisely 


controlling their geometry during microfabrication, we 
have devised technique that lithographically controls 
higher-order dispersion over an octave of bandwidth. 
After introducing the approach and microfabrication 
method, resonator dispersion is measured at various 
bands from 960 nm to 2100 nm and compared with 
modeling to confirm dispersion control. 

Results 

Dispersion design principle. Optical fibre designers 
use multiple cladding layers to control both the magni¬ 
tude, sign and spectral profile of the combined material 
and waveguide dispersiorP^HUl. We devised a method for 
dispersion control in disk resonators that is analogous 
to these methods. Fig. la compares the dispersion of 
a single-wedge disk with that of a single-core optical 
fibre. Details on fabrication and basic properties of the 
single-wedge disk are contained in ref [43]. Briefly, the 
device is fabricated from thermal silica on silicon using 
lithography combined with wet and dry etching. In 
terms of dispersion control a key parameter is the wedge 
angle 0 . The group index spectra of three resonators 
that differ only in their wedge angles are shown in the 
figure. The group index of bulk silica is also provided 
for comparison. In effect, the wedge angle introduces 
a component of norma l dis persion that is stronger 
for smaller wedge angleBy comparison, in the 
single-clad fibr^, the group index is varied (for fixed 
core diameter) through control of An. 

Fig. lb shows the calculated optical mode in the 
single-wedge disk and the optical fibre to further un¬ 
derstand the mechanism of dispersion control in each 
system. In the resonator, longer wavelength modes 
see their centroid of motion around the resonator 
shift inwards (i.e., smaller radii) so that the effective 
optical p ath is smaller. This is equivalent to normal 
dispersioiP^^^. In the fibre, longer wavelength modes 
have a greater spatial overlap with the lower-index 
cladding. Greater contrast between the core and 
cladding indices therefore strengthens this normal 
contribution to dispersion. 

This similarity between the resonator and fibre cases 
suggests a method to engineer dispersion in the resonator 
that is illustrated in the Fig. Ic. In the right-hand 
panel, the case of a multi-clad fibre is considered in 
relation to a resonator featuring multiple wedge angles. 
The multi-clad fibre can be understood by analyzing the 
dispersion in short and long-wavelength case^. Short 
wavelength modes are confined primarily by the core 
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FIG. 1: Fibre-inspired cavity dispersion design, (a) Upper panels depict single-wedge resonator and single¬ 
core fibre; and lower panels give the corresponding group-index spectra of the fundamental modes in the single¬ 
wedge disk (left panel) and single-clad fibre (right panel). Here, 0 is the wedge angle of the resonator, and An is 
core-cladding refractive index difference. The dashed curves give the bulk group index of silica. Calculations assume 
a resonator diameter of 3 mm and fibre core diameter of 11 jim. Micrograph of resonator cross section and fibre 
refractive index profile are provided as insets. White scale bar is 10 fim. (b) Finite element simulation of the fun¬ 
damental mode in a single-wedge disk (left panel) and single-clad fibre (right panel) at the wavelengths 1000, 1500, 
and 2000 nm. The index contrast and core diameter in the fibre mode calculation has been set to An = 0.0028 and 
4.5 /im so as to make the variation in mode profile with wavelength more readily observable, (c) As in panel a, 
but for double-wedge disk and multi-clad fibre. The dashed lines in the group index spectra give either the single 
wedge results (10° and 45° cases) or the single core and clad cases for the fibre. The blue curve in the lower right 
panel gives the multi-clad fibre case. The blue and red curves in the lower left panel give the multi-wedge resonator 
in cases where the angle ordering is red: 45° (outer)->10° (inner) and blue: 10° (outer)^45° (inner), (d) Same as 
panel b, but for a double-wedge resonator and multi-clad fibre. The core diameter in the multi-clad fibre plots are 
set to 4.5 /rm to make the variation in mode profile with wavelength more readily observable. 


region and experience a group index that is similar to 
a corresponding single core fibre. Longer wavelength 
modes are primarily confined by the secondary cladding 
layer. The corresponding spectral dependence of the 
group index can be understood as a transition between 
these two extreme cases (see spectrum for multi-clad 


fibre in Fig. Ic). The wavelength at which the tran sition 
occurs can be correspondingly engineere d^^ F^ For 
example, by placing the outer, higher index region closer 
to the inner core, the transition will occur at shorter 
wavelengths. 

Using analogous reasoning for the resonator, the case 
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FIG. 2: Microfabrication process flow and side-view micrographs, (a) Fabrication process flow steps re¬ 
quired to fabricate a double-wedge resonator, (b) Measured wedge height (“t”, see inset of d) control as a function 
of etch time, (c) Measured Q factor of TM fundamental mode in single, double, and quadruple wedge disks, (d) 

D 2 (=AFSR per mode) as a function of the outer wedge height (“t”, see inset) of double-wedge disk, (e) as a 
function of the third wedge height (“d”, see inset) of quadruple-wedge disk. Outermost wedge height is designed to 
make D 2 within 0-1 kHz (see upper panel), (f-h) Scanning electron micrographs of side-views for (f) single, (g) 
double and (h) quadruple wedge resonators. 


of a double-wedge resonator is considered in the left-hand 
panel of Fig. Ic. In this case, the shorter wavelength 
modes experience the outer, larger-angle wedge (^i=45°) 
while the longer wavelength modes experience the inner, 
smaller wedge angle (6^2=10°). The spectral dependence 
of the group index thereby transitions between these 
two extreme cases. By controlling the specific wedge 
angles the long and short wavelength limiting behavior 
is controlled. Likewise, in analogy with the multi-clad 
fibre, the location of the wedge-angle transition allows 
control of the wavelength band over which the group 
index transitions between these two limits. As shown 
below, this concept can be extended to the design the 
group index of a quadruple-wedge disk which provides 
more flexibility in dispersion control over a wider range 


of wavelengths. 

Device structure and fabrication process. The 

multi-wedge disk process is an extension of the tech¬ 
niques developed for the original single-wedge disk 
resonatoJ^. Fabrication begins with thermally grown 
Si02 on Si, which is then processed with photolithogra¬ 
phy and wet etching into single-wedge shaped disks (see 
top row in Fig. 2a) in which the wedge angle can be 
controlled from 8° to 55° (see Methods). A double-wedge 
(see lower row in Fig. 2a) is created through another 
cycle of thermal oxidation, photolithography and wet 
etching. During this second cycle, the angle of the 
additional wedge is adjusted as before and its height is 
controlled through the duration of etching. Calibrations 








































4 


RF source (f ) 

' eo' 



resonator 


b 



resonator 


DBR 
Laser 
(895 nm) 


NL 

LPF 


FIG. 3: Dispersion measurement setup, (a) Broadband (1400—1700 nm) FSR measurement setup using ex¬ 
ternally broadened EOM comb. ECDL: External cavity diode laser, EDEA: Erbium doped fibre amplifier, HNLE: 
Highly nonlinear fibre, TE: Tunable bandpass filter, PD: Photodetector, (b) ESR measurement setup at 2100 nm 
using difference frequency generation between a 895 nm distributed Bragg reflector (DBR) laser and 1550 nm EOM 
comb. NL: Nonlinear crystal, LPE: Lowpass optical filter 


have shown that the height increases at a rate of ^ 10 
nm/min (cf. Eig. 2b), resulting in the following spectral 
rates: ESR ^ 400 kHz/min and AESR ^ 300 Hz/min 
(cf. Eig. 2d). The process can be repeated to add 
additional wedges. A four wedge structure is shown in 
Eig. 2 h. Because of the added complexity of the process, 
the optical Q factors of the multi-wedge structures are 
lower than for the single-wedge devices. Nonetheless, Q 
factors of 1.63, 1.12, 1.62 x 10^ for single-, double-, and 
quadruple-wedge, respectively were achieved at 1550 
nm. 


Dispersion characterization. To characterize dis¬ 
persion in the resonator we measure FSR^ {fi is 
mode index) in several wavelength bands and calculate 
^ 2 (m) = dFSR ^'/which is related to the GVD 
parameter through the following expression: 


/32(m) 



1 D2{^J) 

' FSm 

fj, 


( 1 ) 


where R is a resonator radiui^^^. As a preliminary test 
of D 2 (at 1550 nm) control in double-wedge design, a 
series of the outer wedge heights (“t”, see inset) were 
fabricated with different wet etch times. A plot of the 
measured D 2 versus outer wedge height is provided in 
Eig. 2d, and the solid line is calculated using a finite ele¬ 
ment solver. D 2 increases as a function of etch duration 
at the average rate of 300 Hz/ min, and the measured 
D 2 has a standard deviation of 495 Hz from the finite- 
element-simulation resulting from the available precision 


in the atomic force microscope and the scanning electron 
microscope measurement of the resonator cross-section. 

An additional interior wedge adds geometric disper¬ 
sion (normal dispersion) in the spectral range that opti¬ 
cal modes experience the new wedge angle. Eig. 2e shows 
measurements of D 2 and D 3 in such a structure plotted 
versus the parameter “d” (see inset in Eig. 2e). Here, we 
define higher order dispersion term as 

//X = /o + • /^ + - 

where /o is the frequency at which the dispersion is 
measured and Di, D 2 , D 3 correspond to ESR, AESR 
per mode, and the third order dispersion paramete j^^ * ^^ l 
As shown in measurements below, the slope of dispersion 
curve at 1550 nm gets flatter (D 3 « 0) as “d” decreases 
and eventually becomes inverse (D 3 > 0 ) to the slope of 
material dispersion curve. In addition, the outer wedge 
height is also designed to keep D 2 consistent within a 
range of 0 — 1 kHz. 

To measure the ESR quickly and over a broad range 
of wavelengths, we have modified a method reported 
in ref [52]. The measurement setup is shown in Eig. 
3a and incorporates a frequency comb generated using 
electro-optical modulation (EOM) in conjunction with 
optical broadening using high nonlinearity fibre (HNLE). 
The seed laser for the EOM comb was a tunable laser 
that was scanned continuously during the measurement. 
Eurther details on this EOM comb are provided in 
ref [53|. The comb bandwidth was approximately 300 
nm centered around 1550 nm and the line spacing 
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FIG. 4: Dispersion measurements for single-, 
double-, and quadruple-wedge disks. FSR, and 
D 2 for (a) single-, (b) double-, (c) quadruple-wedge 
disk. Blue and green data points are measured FSR 
and D 2 , respectively; red and orange lines are numeri¬ 
cally calculated FSR and D 2 , respectively. The cavity 
geometry is directly imported from the scanning elec¬ 
tron micrograph image (inset, bar: 10 jim) for finite 
element simulations. 


was tunable and set to coincide approximately with 
the resonator FSR. A relatively narrow band of comb 
frequencies was coupled to the resonator using a tunable 
optical filter and the transmitted optical signals were 
measured using a photodetector. A measurement 
proceeded by first setting the tunable filter to a desired 
wavelength. The optical filter had a 3-dB-bandwidth of 
3 nm so approximately 15 comb teeth would be launched 
to the resonator input port (the EOM comb line spacing 


was approximately 22 GHz). The photo detected signal 
was then observed as a time trace on an oscilloscope 
that was synchronized with the scanning seed laser. If 
the microwave frequency used to establish the EOM 
comb line spacing was set to exactly coincide with the 
ESR of the resonator at the selected wavelength, then 
a single transmission minimum would appear since the 
filtered comb teeth are all coupled to the neighboring 
resonances of the resonator. By slightly tuning the 
comb drive frequency away from the ESR value, the 
single transmission minimum would break up into many 
separate peaks since the comb teeth would now achieve 
resonance at slightly different points in the laser scan. 
The resolution in determination of the ESR by this 
method was estimated to be about 100 kHz. It is also im¬ 
portant to note that the method provides an ESR value 
that is spectrally averaged since multiple comb teeth 
participate in the process. This fact tends to smooth 
the data with respect to effects like avoided-crossings of 
mode families. In addition to measurements in the 1500 
nm band, the EOM comb was also shifted in wavelength 
to 2.1 /im by difference frequency generation using a 895 
nm pump laser. Also, measurements of ESR using the 
technique reported in ref [52] were conducted at 980 nm 
and 1300 nm. Collectively, these measurements provided 
a very good picture of how dispersion could be modified 
through judicious cavity design. 

Eig. 4a-c show measured ESR, D 2 , and the results 
of simulations for the device structures provided in the 
insets. As has been noted in ref [52|, a single wedge 
contributes geometrical dispersion that is both normal 
and approximately spectrally flat. Also, smaller wedge 
angles increase the strength of this contribution. In 
figure 4a, the zero dispersion wavelengt h (Xz pw)^ has 
been shifted by this effect to 1.5 (material 

dispersion has a Xzdw near 1.3 iim). 

In the double wedge resonator (Eig. 4b) there is 
a gradual onset of the geometrical dispersion from 
shorter to longer wavelengths as a result of a steep outer 
wedge (low normal dispersion) to the shallow interior 
wedge (larger normal dispersion). The overall effect is 
to provide a flattening of the dependence of ESR on 
wavelength over a broad range of wavelengths. The 
dispersion in the transition is therefore reduced. The 
addition of additional wedges can be used to extend 
control to longer wavelengths. In the quadruple-wedge- 
design the additional interior wedges provide control out 
to 2 fim (Eig. 4c). 

Conclusion 

This study has proposed and tested a resonator dis¬ 
persion control approach that is inspired by dispersion 
control methods in multi-cladding optical fibre. This ap¬ 
proach provides the first method of controlling resonator 
dispersion over a broad range of wavelengths (equivalent 
to having control over high-order dispersion). The 
approach builds on the high-Q silica-on-silicon wedge- 
resonator fabrication by adding additional wedge angles 
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to the resonator. The position of each wedge and its 
angle determine the wavelength location and strength 
of a corresponding contribution to normal dispersion. 
As a result, the approach is highly intuitive and when 
combined with numerical modeling provides a well con¬ 
trolled way to engineer dispersion over broad wavelength 
bands. Dispersion engineered high-Q resonators, as 
described here, will be useful in nonlinear resonator 
devices such as frequency microcombs. Moreover, the 
methods can also be applied in control of dispersion in 
nonlinear waveguide structures. 

Methods 

Device fabrication. The fabrication process was based 
on the previously described method for the chemically 
etched wedge-resonatoiE^. Disk resonators were fab¬ 
ricated on (100) prime-grade float-zone silicon wafers. 
The first Si02 layer was thermally grown at 1000° C with 
a thickness in the range of 7—9 /im, and photoresist was 
patterned using a KarlSuss MA-6 aligner on the oxide 
layer. Hexamethyldisilazane (HMDS) was additionally 
deposited prior to photoresist coating for wedge angles 
larger than 20° (right column of Fig. 5a). For wedge 
angles in a range of 8°—20°, photoresist was spin-coated 
on the oxide layer without HMDS (left column of Fig. 
5a). Photoresist thickness was used to adjust the wedge 
angle. The photoresist pattern acted as an etch mask 
during immersion in a buffered hydrofluoric acid (HF) 
solution, and the resist was removed once wet etching 
was completed. The second Si02 layer was grown at 
the same temperature as the etched wedge disk, and 
the thickness of the oxide was accurately controlled via 
duration of thermal oxidation. Again, the photoresist 
was patterned on the oxide layer and the wet etch 
proceeded to create additional wedge angles. After 
cleaning to remove photoresist and residues, XeF 2 dry 
etch was applied to etch the Si. 

Finite element simulation. The numerical calcula¬ 
tions were implemented using the COMSOL multiphysics 
package. Fibre parameters used in the Fig.l were from 
ref iHnasi, and the values for fused silica refractive 
index were taken from the Sellmeir equation of ref ESI- 
Resonator dispersion was calculated via an iterative 
approaclP^. Fig. la-b use a resonator thickness of 8 
/im, and Fig. Ic-d use outer wedge heights of 3 i^m and 
5 /rm for the angle ordering 45° —10° (red) and 10°—45° 
(blue), respectively. Fig. 4 uses a resonator geometry 
obtained from fabricated devices. Oxide thickness was 
measured using Filmetrics model F40 and Dektak 3ST 
profilometer, and the wedge profile was obtained using 
an atomic force microscope (Bruker Dimension ICON) 
and scanning electron microscope (Hitachi S-4100). Fig. 
5b shows an atomic force microscope image of top and 
bottom surfaces of double-wedge and quadruple-wedge 
disks. The top surface profile was obtained from the disk 
before silicon dry etch (5th step in Fig. 2a), and then 
the silica disk was completely removed by BHF and the 



FIG. 5: Wedge angle control and atomic force 
microscope image of double and quadruple 
wedge structures, (a) Scanning electron micro¬ 
graphs of cross sections for single wedges. The left col¬ 
umn {0 =8, 9.5, 11, 14°) shows cross sections of oxide 
wedges processed without HMDS, and the right column 
{0 =18.5, 23.5, 27, 55°) shows the structure processed 
with HMDS. Within each of these cases, the wedge an¬ 
gle can be accurately controlled by adjustment of pho¬ 
toresist thickness, (b) Atomic force microscope image 
of double wedge and quadruple wedge (top and bottom 
surfaces are superimposed). 


silicon surface morphology was scanned for the bottom 
surface profile. 

Dispersion parameters fitting. The resonance 
frequencies of one mode family can be described as a 
Taylor series shown in eqn.(2), and the coefficients of 
the series correspond to dispersion parameters at /i = 0. 
The FSR can also be expressed as a Taylor series with 
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FIG. 6: Dispersion parameters fitting (a) Mea¬ 
sured FSR of single-wedge disk and the Taylor series 
fit (solid line). Blue dots represent the measurement 
value using the EOM comb method (Fig. 3 a), an d red 
dots are measured by sideband spectroscop}^^^^ for 
comparison. The FSR is 22.042 GHz at 1550 nm, and 
dispersion parameters are determined as the coefficients 
of Taylor series at the same wavelength, (b-c) Same 
as panel a, but for quadruple-wedge disks. FSR: 21.36 
GHz for (b), 21.40 GHz for (c) at 1550 nm. 


the same coefficients in eqn.(2). 

FSRifi) = + £>2 ■ M + + ■ ■ ■ (3) 

Fig. 6 shows the measured FSR of single-wedge and 
quadruple wedge structures (Fig.4), and the fitted 
polynomial curves (eqn.(3)) to the experimental results. 
Here, FSR was measured using both EOM-assisted 
(Eig. 3a, 300-nm-bandwidth) and sideband^ (70-nm- 
bandwidth) spectroscopy methods. The variation of 
ESR is very small (<10 kHz /mode), thus it is neces¬ 
sary to accumulate the variation of ESR throughout 
multiple fseP in order to fulfill the condition that 
I D 2 I -/i > A/, I ^ I > A/ where A/ is the 
resolution in determination of the ESR (approximately 
100 kHz). 
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